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1. Introduction. 



There are many phenomena in nature which can be considered as some manifestation 
of blow ups, like hurricanes, tornadoes, sandstorms, etc. If we believe that Navier-Stokes 
system describes well enough the motions of real gases and fluids under normal conditions, 
then it gives some reasons to expect that blow ups in solutions of this system also exist. 

We consider in this paper the 3£>-Navier-Stokes system for incompressible fluids moving 
without external forcing on R 3 with viscosity equal to 1. After Fourier transform it becomes 
a non-local, non-linear equation for a non-known function v(k,t) with values in C 3 , k G -R 3 , 
t > 0. The incompressibility condition takes the form* (v(k,t), k) = and 

t 

v(k,t) = exp{-t\k\ 2 }v(k,0) + i J exp{-(t - s)\k\ 2 } ds- 

(1) 

/ (v(k-k',s), k) ■ P k v(k',s)dk' 

In this expression v (k, 0) is an initial condition and P k is the orthogonal projection to the 
subspace orthogonal to k, i.e. P k v = v— ^°/^ k \ k ■ The formula (1) shows that the Navier-Stokes 
system is genuinely infinite-dimensional dynamical system: the value v(k,t) is determined 
by the integration over all "degrees of freedom" and previous moments of time. 

The problem of blow ups in solutions of the Navier-Stokes System(NSS) appeared after 
classical works of J. Leray (see [Le 1]) where he proved the existence of the weak solutions 
of NSS. O. Ladyzenskaya proved the existence of strong solutions of 2-dim NSS in bounded 
domains (see [La 1]). Many important contributions to the modern understanding of the 2- 
dim fluid dynamics were done by C. Foias and R. Temam (see [FT]), V. Yudovich (see [Yl]), 
Giga ([Gl]) and others. However, the situation with the 3-dim NSS remained unclear. The 
Clay mathematical institute announced the problem of existence of strong solutions of the 
3-dimensional NSS as one of the most important problem in mathematics of the XXI-century 
(see [CI]). 

In this paper we omit the condition that v(k,t) is the Fourier transform of a real-valued 
vector field in the x-space and consider (1) in the space of all possible complex-valued func- 
tions with values in C 3 . In this situation the energy inequality does not hold. Detailed 

* Since k G R 3 , v(k,t) e C 3 , the order in the inner product is important. 



assumptions concerning the initial condition v(k, 0) will be discussed later (see §7). In all 
cases v(k,0) will be bounded functions whose support is a neighborhood of some point 
(0, 0, The incompressibility condition implies that the components v i(fc, 0),v 2 (k, 0) of 

v(k, 0) are arbitrary functions of A; while t> 3 (A;,0) can be found from the incompressibility 
condition (v, k) = 0. 

Various methods (see, for example, [K], [C], [SI]) allow to prove in such cases the existence 
and uniqueness of classical solutions of (1) on finite intervals of time. For these solutions 
(see, for example [S2]) 

\v(k, t)\ < const exp {— const Vt ■ \k\} , < t < t . (2) 

Presumably, v(k,t) has an asymptotics of this type but this requires more work. According 
to a conventional wisdom, possible blow ups are connected with the violation of (2). 

In this paper we fix t and consider one-parameter families of initial conditions i>a(/c, t) = 
Av(k, 0). We show that for some special v(k, 0) one can find critical values A cr = A cr (t) such 
that the solution VA cr (k, s) blows up at t so that for t' < t both the energy and the enstrophy 
are finite while at t' — t they both become infinite. Even more, for t' < t the solution 
decays exponentially outside some region depending on t. As t' f t this region expands to 
an unbounded domain in R 3 . 

Our main approach is based on the renormalization group method which is so useful 
in probability theory, statistical physics and the theory of dynamical systems. It is rather 
difficult to give the exact formulation of our result in the introduction because it uses some 
notions, parameters, etc., which will appear in the later sections. Loosely speaking, we show 
that in ^-parameter families of initial conditions, for £ = 10, one can find values of parameters 
for which the solutions develop blow ups of the type we already described. The meaning of 
£ is explained in §4, §5, §6. 

We thank C. Fefferman, W.E, K. Khanin and V. Yakhot for many useful discussions. A 
big part of the text was prepared during the visit of the second author of the Mathemat- 
ics Department of California Institute of Technology and we thank the Department for its 
very warm hospitality. We also thank G. Pecht for her excellent typing of the manuscript. 
The financial support from NSF Grant DMS 0600996 given to the second author is highly 
acknowledged. 
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2. Power Series for Solutions of the 3"D-Navier-Stokes-Systems and 
Preliminary Changes of Variables 



Our general approach is based upon the method of power series which were introduced 
in [SI], [S2]. We write down the solution of (1) in the form: 

t 

v A (k,t) = exp{-t|A;| 2 } • Av(k,0) + / exp {-(t - s)\k\ 2 } • A p h p {k, s) ds (3) 



The substitution of (3) into (1) gives the system of recurrent equations connecting the func- 
tions h p : 

h!(k,s) = exp{-s|A;| 2 }t;(A;,0), (4) 
h 2 (k,s) = i J (v(k-k',0),k)P k v(k',0) ■ exp {-s\k - k'\ 2 - s\k'\ 2 } d 3 k' , (5) 

R 3 

h p (k,s) = ij ds 2 J (v(k - k',0), k) Pkh p -i(k' , s 2 ) ■ 

PS PS P 

exp{-s\k - k'\ 2 - (s - s 2 )\k'\ 2 } d 3 k' + i ^ ds 1 ds 2 (h Pl (k - k' , si), k) ■ 

„,+ m =„Jo Jo J 



P1+P2=P 3 
P1>P2> 1 n 



Pkh P2 (k', s 2 ) • exp{-(s - Sl )|A; - kf - (s - s 2 )\k'\ 2 } d 3 k' + 

l So dSl I ^ ^"^'^'^ Pkv(k',0) ■ exp{-(s- Sl )\k- k'\ 2 -s\k'\ 2 }d 3 k' . (6) 



R3 

Clearly, h p (k, s) _L k for every p > 1, k e R 3 . 

It follows from the results of [S2] that the series (3) converges for sufficiently small s and 
gives a classical solution of (1). Make the following change of variables which simplifies (4), 
(5), (6). Put k = ky/s, k! = k'y/s, introduce relative times §i,s 2 ,Si = §is, s 2 = s 2 s and 
denote g r (k, s) = h r , r > 1. Then 

9l (k,s) =exp{-|fc| 2 }- v(±0] , (4') 



g 2 (k,s) = h 2 | A s j = A. J ( v f ) 



R 3 



P s u |^,0] exp {-|jfe - fcf - \~k'\ 2 } d 3 ~k' , (5') 



i 



g P (h,s) = l - j ds 2 J (v I k fc , j , k ) • P- k 9 P ^i(ky^s 2 ',S2s) 
exp{-|jfc-fc , | 2 - (1 — S 2 ) l^'l 2 }^ 3 ^' + 



l l 



+ i ^ y d§i / ds 2 / (g Pl ((k - k') \Zh,S!s),k) 

Tl, = J J J 



P1 >l,p 2 >l U U R 



P~ k 9p2(k'Vs2, hs) exp {-(1 - h) \k - k'\ 2 - (1 - s 2 )\k'\ 2 } d 3 k' 



~s J dSl J ^ p_1 ^ fc ~ ^)) SlS ' ) ' ( 



+ 

s 

R3 



exp{-(l-5 1 )|A;-A; , | 2 - \k'\ }d k' (6') 

The function g 2 (k, s) has a singularity at s = even in the case of functions with compact 
support: for small s its values are of order —=. This singularity is integrable and all g p (k, s), 

p > 2, are bounded. The singularity is connected with our choice of the coordinates k, k'. 

The formulas (4)-(6) or (4')-(6') resemble convolutions in probability theory. For example, 

if C — supp-u (k, 0) then supp h p — C + C H h C. Therefore it is natural to expect that 

p times 
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hp and g p satisfy some form of the limit theorem of probability theory. This question will be 
discussed in more detail in the next sections. 



Make another change of variables. Assume that we have some p. The terms in (6') with 
Pi < P 1 ^ 2 and P2 < p 1 / 2 will be called boundary terms. They will be treated as remainder 
terms and will be estimated later. Suppose that we have some number which later will 
be assumed to be sufficiently large. Introduce the vector K,^ = (0, 0, rk^). These will be the 
points near which all g r will be concentrated, 

p 1 / 2 < r < p — p 1 / 2 . We write k = K.^ + \fr • Y, Y e R 3 . Thus instead of k we have the 
new variable Y = (Y 1 ,Y 2 ,Y 3 ) which typically will take values 0(1). Put «(°) = (0,0,A; (0) ). 

In all integrals over s±, §2 in (6') make another change of variables 1 — s j — jp, j = 1, 2. 
Instead of the variable of integration k' introduce Y' where k' = + y/pY'- We write 
g r (Y, s) = g r (KV + ^Fy, s), 7 = f, f = 1 - 7- Then from (6') 



g p (Y,s) = g p (K» + JpY,s) = p*' 2 



pi vl 



P1,P2>Vp 
P1+P2=P 







9pi 



R3 



pi ■ P2 

Y' 



exp < —Vi 



Y-Y' 



Y' 



v^(!-7) 



to) + y g P2 , r 



d 3 Y' 



(7) 



This is the main recurrent relation which we shall study in the next sections. It is of some 
importance that in front of (7) we have the factor p 5 ^ 2 and inside the sum the factor 4^ • 4?. 
Both are connected with the new scaling inherent to the Navier-Stokes system. 



j3. The Renormalization Group Equation 



As p — > 00 the recurrent equation (7) takes some limiting form which will be derived in 
this section. All remainders which appear in this way are listed and estimated in §8. 



The main contribution to (7) comes from pi,p 2 of order p. If Y, Y 1 



0(1) then 



are small compared to = (0, 0, A;^ ^). Therefore the Gaussian term in (7) can be replaced 
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by exp{ — (9 1 + 9 2 )\k^\ 2 }, si and s 2 can be replaced by 1 and the integrations over 9±, 9 2 
and Y' can be done separately. Thus instead of (7) we get a simpler recurrent relation: 



Pl>P2>P 1 ^ 2 * 



•P s( o, + y.fe I rf 3 r. (8) 



v / a 3 7)'' 



In view of incompressibility 



= £ (<7 P1 (^ {0) ' Pi + ^ ' v^ s ) > + ^ " 7v^ 

+ if <fe(^,s),*V(l^)> (9) 



Write p p in the form 



4 P (E, s) = (G«(y, s) , Gf } (Y, s) , -L fW(y, s )) . (10) 



Since fc = ■ p + Yy/p, the incompressibility implies 



~ k 
< g r (k,s), k > = < g r (k,s), - >= (11) 



and for Y = 0(1) 

4 • ^CM + ^=G< r) (r,*) + ^ • F«(F )S ) = O (- 

y/r Vr ^ \r 



(12) 



In our approximation we replace (12) by 



Y X (#\Y,8) + Y 2 GP(Y,s) + F^(Y,s) = 0. 



(13) 



Thus for given Y 1: Y 2 , Y 3 the component F r can be expressed through G[ , G 2 . This remains 
to be true even if we do not neglect the rhs of (13). Return back to (9). From (13) 



7-1 

V7 



< 9 P 



1 { Y J ,S ) ' 



Y-Y' 



> 



i 



7_ 



+ 



y 2 ' (^(pO / y-y 



v^a v 7 ! 3 ^ V Vy J J 



; l-7 2 V v 7 ^ 

In our approximation the inner product in (14) can be replaced by 

7 V v 7 ^ 1 \ \fi ) 2 \ V7 ' / / 

/ v . . /v v \ 



l 



V^7 



' 



VI -7 2 



v/7 



According to the definition of the projector 
P«(o )+ ^fe(^,s) = 9 P2 (^=, s ) 



(14) 



(15) 



< 9 P 2 



<^(o) + y « ( o) + y > 



9p2 



Y' 



+ o 



(16) 



This shows that in the main order of magnitude the projector is the identity operator and 
we come to a simpler recurrent relation instead of (8): 
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g p (Y,s) = i\i ■ 



p z 



P1.P2 

>pl/2 
P1+P2=P 



pi ■ vl Jr? 



7-1 fY 1 -Y{ 



(pi) 



Y-Y' 
~7T 



Y2 ~ Yi g(pi) 



a/7 



Y-Y' 



a/7 



,s + Vl^7 



*7 ^) ^ - Y 



a/7 



,s + 



V 

I 2 q(pi) 



Y-Y 1 
"TP" 



• #P2 



y 



(17) 



The main assumption which we shall check below in the next sections concerns the asymptotic 
form of g p (Y, s) as p — > 00: for some interval S&> = [S {p) , S {p) ] on the time axis and some 
A, positive and for all r < p 



g r (Y,s) = A'-V 



a 



(i) 



a 



(i) 



(2) 



( . xp i — (\ Yl \ 2 + \Y 2 2 \) l ■ yi— exp 
(H 1 (Y 1 ,Y 2 ,Y 3 ) + S[ r) (Y,s), H 2 (Y u Y 2 ,Y 3 ) + S 2 \y, s) , S 3 \y, s)) 



a 



(2) 



(18) 



where <5j r) (Y» — ► as r — > 00, j = 1, 2, 3. Later we shall explain in more detail in what 
sense the convergence to zero takes place. The substitution of (18) into (17) gives 



9 P (Y, s) 



p ■ A p - 



Ei 1 1 
_. 72(1-7)2 



pi 

7= — 

' P 



i? 3 



7-1 



\/7 V >/7 



a/7 



+ 



a/7 

y 



2 #2 



v/7 

y-y 

V7 



a/T^7 



v/7 



if 



y 



27T7 



cxp 





/iy-y 


2 + 


Y 2 - Y 2 ' 




V 7 



r(l) 



2tt(1 -7) 



exp < — 



<7« \y;\ 2 + |y 2 '| 



a 



(2) 



27T7 



cxp 



2 I-7 

^ (2) |y - ^ 



7 



a 



(2) 



2?r(l -7) 



cxp 



(7 



(2) ly/12 



2 1-7 



d 3 Y' 



(19) 



Here 



ff 



y 



#1 



y 
1 2 



y 



Vi — 7' Vi = 7' Vi = 7/ 



ff, 



Vi - 7' Vi _ 7' V 1 - 7 



,0 



We do not mention explicitly the dependence of ff on s. 



The last sum looks like a Riemannian integral sum whose limit takes the form as p — ► 00: 

r(l) 



Aexpl-^lyp + lynj.^.expj- 



a^\Y 



2 ! V 2tt 



H{Y) 



a 



\kW\-.iu 

(i) 



72 (1 — 7) 2 d'j 



a 



(i) 



R3 



27T7 



cxp 



2tt(1 -7) 



exp 



+ \Y 2 f)\ a 



a 



(2) 



2tt(1 -7) 



cxp 



2(1-7) 



gWjYtf 
2(1-7) 



aW(\ Yl -Y(\ 2 + \Y 2 -Yj\ 2 ) 
2 7 

*<a\Y 3 -n\* 



r(2) 



27T7 



cxp 



27 



V7 V V7 



v/7 



y - n 



#9 



y-y r 

V7 , 



+ 72(1 - 7 ) 



y 

v/1^7 



#1 



y-y 
^7" 



+ 



y 



V 1 -7 VV 1 -7 



y-y 



if 



y 
v / r^7 



rf 3 y 



(20) 



The integral over Y 3 is the usual convolution. Therefore we can look for functions ffi, ff2 
depending only on Y U Y 2 , i.e. H^Y) = H 1 (Y 1 ,Y 2 ), H 2 (Y) = H 2 (Y!,Y 2 ). Write down the 
equation for H 1 , H 2 which does not contain Y 3 : 



9 



exp 



Y\ 



a 



(i) 



2tt 



H(Y) 



R 2 



a 



(i) 



27T7 



cxp 



a 



(i)|y _ y 



/|2 



27 



(7 



(1) 



2tt(1 -7) 



exp 



2(1-7) 



7 )S 



3/2 / n-y/ 



Y-Y' 



+ 



y-y 



0) 



(21) 

Here F = (Y u Y 2 ), Y' = (Y{, Yf), H(Y) = {Hx{Y x , Y 2 ), H 2 (Y h Y 2 ). This is our main equation 
for the fixed point of the renormalization group which we shall analyze in the next section 
(see also §7). 

§4. Analysis of the Equation (21) 

The solutions to the equation (21) have a natural scaling with respect to the parameter 
a = Namely, if we solve the equation (21) for a = 1 and let the corresponding solution 
be H(Y), then the general solution for arbitrary o is given by the formula 

H a (Y) = ^H(^Y). (22) 

This is analogous to the usual scaling of the Gaussian fixed point in probability theory. 
Thus, it is enough to consider the equation (21) for a — 1. We shall show that there exists a 
three-parameter family of solutions to the equation (21) for a — 1. The equation (21) takes a 
simpler form if we use expansions over Hermite polynomials. All necessary facts about Her- 
mite polynomials are collected in the Appendix 1. For H{Y 1 ,Y 2 ) = (Hi{Yi,Y 2 ), H 2 (Y 1 ,Y 2 )), 
we write 

H j (Y 1 ,Y 2 )= h^(m 1 ,m 2 )He mi (Y 1 )He m2 (Y 2 ), j = 1,2 (23) 

mi,m2>0 

where He m (z) are the Hermite polynomials of degree m with respect to the Gaussian density 
^= exp {-\z 2 }. We have (see If33jl ): 

zHe m (z) = He m+ i(z) + mi/e m _i(z) , m > (24) 
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and 

Heo(z) = 1, zHeo(z) = z = He\(z) 



Also we use the formula (see 

:y-y'\ i 

He r 



7 



\Y - Y'\ 2 
exp <; — ^> //• 



111-2 



V 



\A -7/ v 7 ^ 



exp 



/|2 



2(1- 7 )J ^ =7 2 ^~ 7 ) 2 He m 1+ m2(X) ^J= exp j } 



(25) 



Substituting (23) into (21) and using (24), (25), we come to the system of equations for the 
coefficients h(mi,m 2 ) which is equivalent to (21): 

h^( mi ,m 2 ) = J2 J $m» ■ {(5ifc (1) )(miX) + (« (2) )KX 2 )} fe^K.m*) 



11 " 



+J2m» ■ {h {1) {m' 1: m' 2 ) {B 1 h^){mlmD + m' 2 ) (^K.rnJ) } (26) 



where m' = + m 2 , m" = m" + m' 2 ' and 



J (1) 

J m'm" 



t(2) 
^m'm" 



7 2(1—7) 2 (£y 



m_+l /, s rn"+2 

72 (1 — 7) 2 dry 



(27) 



(B\h^) (m\,m 2 ) = h^' (m^ — 1, m 2 ) + (m' x + l)/i^'(m^ + l,m 2 ) 



(5 2 /i (j) ) (morn's) = /i^ {m' l ,m 2 - 1) + (m' 2 + 1) h u) (m' 1 ,mf 2 + 1) 



To simplify the system (26), we shall look for solutions with h^\0,0) = 0, j = 1,2. 
Below we sometimes write (mi,m 2 ) as /im| m2 or hm ltm2 when there is no confusion. 
Similar conventions will be applied to J^(m 1 , m 2 ). For mi + m 2 = 1, we have 
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r h (i) 

'no 


- J (1) 

— J 01 




+ 4?) 


'no 


+ J (2) 

+ J 10 


(h (1) h (1) 

{'no 'no 


-t- il 10 n m ) 


h {1) 

"Ol 


- J (1) 

— J 01 




+ 4?) 


'^01 


+ J 10 


\"oi 'no 


+ h (2) h {1) ) 
-t- /i 01 /i 01 j 


< 

h {2) 
'no 


- J (1) 

— J 01 




+ 4?) 


/i (2) 
'no 


' J 10 


(h il] h {2) 
I'no 'no 


+ "SO 


\ "01 


- 7 (1) 

— J 01 




+ 4?) 


' i 01 


T J 10 


I'nn 'no 


+ CC) 



where J, 



g> = -1/3 and J x {2) 



Case 1. h$ + 4 2) 



01 



1/6. There are two cases: 

6. In this case (4o\ 4i\ 4o\ 4?) om y nee ds to satisfy: 
(4 1 ) +3) 2 = 9-4 1 > (2) 
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This is a two parameter family of solutions. 



Case 2. h$ + h$ ^ -6. In this case {h\Q , , , ) can be uniquely determined and 



,(2) 



(1) 1,(1) 1,(2) ^ 



we have /i 



(i) 
10 



(2) 
01 



-2 h (1) 
z > ' i oi 



(2) 
10 



0. 



For the rest of this paper we shall consider only the case 2 for which h 



(i) 
10 



h 



(2) 
01 



-2, 



h-01 = h\o = 0. Let us write down the recurrent relations for m\ + m 2 — 2, j — 1, 2 : 

/ion — — (2^20^ "I" ^^02^ ^^11^)^20^ ^"^11 ' 4o^ ' 4o^ 4o^ ' ^11 ' ^11 



,(2) 



°20 

{ fcW = 



(J) 
02 



(24? + 4Jo? + 4^)4? + J$>h® ■ (2/&> + 4?) + J{? • 4> • (41 ) + 24?) 



(2<^2o' ) 

-(2J 2 ( ?+ u;;u ir.ru>- 



2 7 {1) • h U) ■ h {2) 



h ij) • J (1) • h {1) 
"01 J n "u 



It is not difficult to check that the only solution to the above system is 4o* 



(?) 

02 



= 0. Solving the recurrent relations for m\ 



m 2 = 3 gives us: 



r h (i) 

"03 


- h {2) - 

— n 30 — u 


, "12 


— "03 


h (1) 
"21 


— "12 


h (1) 
K "30 


— "21 



This shows that (4?; 4i\ 4o) can t> e considered as free parameters. For any p > 4, 
the recurrent relations for m\ + m 2 = p form a linear system of equations for the variables 
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{^mi,p-mi}mi=o with coefficients depending on h^} and /i^ only. In principle, they can be 
solved and an explicit expression for the solutions can be found. We emphasize here that if 
the free parameters take real values then the whole solution is also real. 



l zi i ^30 ) > one can fi 11 ^ an h 



00 

mi ,m2 



It is not difficult to check that for any values of (h^ 2 \ h 
(mi + m 2 > 4) by using (26). The solution we obtain is formal in the sense that it satisfies 
(26) but h mi>m2 with m 1 + m 2 = p may not decay as p — ► oo. We are now ready to formulate 
the theorem concerning the existence of formal solutions to (26). 

Theorem 4.1. For any values of (h^ 2 , ^4i\ ^4o )> there exists a unique formal solution to 
the recurrent equation (26). 



Thus, theorem 4.1 claims the existence of a three-parameter family of solutions of (21) 
parameterized by h^ 2 , h 2 } and ■ It turns out that if h^ 2 , h^} and h^j are suffi- 
ciently small, then h m \, m2 decay as m x + m 2 = p tends to infinity. Let us say that 
h m \,m 2 has degree d if mi + m 2 = d. For each d > 4, introduce the vector = 
(h^ d , • • • , hfl, hfl, . . . , hfl) T . The vector h^> contains all terms of degree d. By 

the recurrent relation (26) 

C (d) h (d) = b (d) ( 28 j 
where the vector contains terms of degree < d— 1. Also £ ]g>(2c(+2)x(2c!+2) j g a ma t r j x: 



C 



(d) 



16d- 16 + 32A; 
~ (d+ l)(d + 3)(d + 5)' 

80d + 80 - 32A; 
~ (d+l)(d + 3)(d + 5)' 

32(d- k + 2) 
'(d+l)(d + 3)(d + 5)' 

32(A; — d — 1) 



(d+l)(d + 3)(d + 5)' 



I 0, 



if l < jfe = Z < d + 1 



if d + 2 < A; = l<2d + 2 



X2<k<d+l,l = d + k 



i(d + 2<k<2d+l,l = k-d 



all other cases 



It is easy to check that if d > 4, then is nonsingular and as d — ► oo, converges to 
the identity matrix. This observation immediately implies the following lemma: 
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Lemma 4.2. Let (C^)" 1 be the inverse matrix of for d > 4. There exists an absolute 
constant C\ > such that for all d> 4 

|| (C^)- 1 || < d . 

We are now ready to derive an estimate which gives the decay of solutions of the recurrent 
relation (26). 

Theorem 4.3. // \ h^\ < 5, \h^i\ < S, \h^j\ < 5 and 5 is sufficiently small, then for some 
C 2 > ; < p < \, we have 

„mi+m,2 

\h^ um2 \<C 2T ^^ Vmi > 0, m 2 > 0, j = 1,2. 

Proof. We begin by noting that hm im2 = if mi +m 2 is even. This can be easily proven by 
using the recurrent relation (26) and the fact that h^f = and fom!,m 2 = for mi + m 2 = 2. 
Let < pi < 1, pi will be chosen sufficiently small. We shall use induction on mi + m 2 
where mi + m 2 is odd. According to the induction hypothesis 

m.i+m.2+2 

l*2Ul < r ( mi+r+7 ) 9(mi + m 2 ) (29) 

for every 3 < mi + m 2 < d — 2 where d > L is an odd number and L will be chosen later 
to be sufficiently large. Also g is a function to be specified later. We shall comment on the 
choice of L and verify the induction hypothesis for 3 < m x + m 2 < L later. Let us show that 
the same inequality holds for m x + m 2 = d. Without any loss of generality, let us consider 
j — 1. The case j = 2 is similar. Fix mi and let bm\ be the (mi + l) th component of the 
vector in the equation (28). We now estimate bm\ using the induction hypothesis (29) 
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and the equation (26): 



d-3 



\h(d)\ < >p j(i) 

IV I — 2.^1 | rn',m" 
m'=2 
d— 3 

I m',m" 

m'=4 
(2-2 

+ V lj (2) 
+ \ J m',m" 

m'=3 

(2-4 

+ V Ij (2) 

I m',m" 

=3 

(K 



2 • 



Pi 



rn"+2 
Pi 



m'=3 
+ 12 



(1) 
(2-2 



r(^) r(^) 



(m' + l)#(m' + l)0(m") 



pf +1 



pf +2 



pf+ 2 



pr" +3 



r(zr^) r(^) 



(m' + 1) • g(m' - l)g(m") 



(m' + 1) • (77/ + 1) • g(m')g(m" + 1) 



pf +2 



pf +1 



r(2iz) r( 



J. 



(i) 



J. 



(2) 



d— 2,2 



t(2) 



r(f) 



The last term in the rhs of the above inequality comes from the case where h m ' iTri2 or h m » m » is 
of degree one since the induction hypothesis holds only for 3 < m 1 +m 2 < d — 2. Also in the 
estimation of the first four terms we use the fact that for fixed (m',mi), there are at most 
min{m' + l,m" + 1} tuples of (m', m' 1 ', m' 2 , m' 2 ) such that wl x + m x = m x , m 2 + m 2 = m 2 , 
m' x + m 2 = m' and m' x + m 2 ' = m". By (27), we have 



J (1) 

m'm" 



r(2) 

u m',m" 



r(™i2)r(^) 
r(=*a)r(^) 



r( 



2 / 



and for some constant C3 > 



J, 



(i) 



2,<Z-2 



J (1) 



+ 



J. 



(2) 



d-2,2 



t(2) 
J l,(2-1 



" d 2 
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Therefore 

„, | < ■ > v 2 / ' r ■ ) 2 ( ■ gym' + l)g(m") 

mi I — t-i / a+7 \ / j T-i / m'+8 \ p 7 m"+7 \ y ' 

\ 2 ) m '=2 \ 2 ) \ 2 ) 

2pf 3 ^ r(f)y + i) rgg) 

H , . , ?N ■ > — - — 7 • — ; „' • gym — I) gym ) 

1 ^ 2 / m '=4 L \ 2 I L \ 2 I 

2 P ( +5 ^ r (2^t3) • K + 1) r(^)-(m" + i) „ 

+ f^±zj • E " • r(ir ^8) " • 9(m')g(m" + 1) 

2p f 3 ^ r(^)-K + i) rffg) 

V 2 / m'=3 \ 2 J \ 2 ) 

j>r_ . ft . _ 12 
r(^) <i 2 r(i±s) p? Jl ' 

PI 



-f7feY ' * " ^ • £ 9(m' + l)<?(m") + 9(m' - l)g{m") 

V 2 / \m/=2 m'=4 



J]) <?K)<?K + 1) + £ 9(m')g(m" - 1) + • ^5 • 

m'=3 m'=3 / 1 ° ' Pl 

where C 4 , C 5 are some constants. Now we specify the choice of the function g. Let gym) be 
such that gi = a and 

m— 1 

g(m) = g(p)g(m — p) for m > 1 
P =i 

By the method of formal power series it is not difficult to show that 

g(m) = I . < 2 '" - w . (2aT 

yK J 2 m! v J 
Clearly, we have const < 9<y ^^ < const, and this immediately gives us 

d+2 ^ d-j-2 

K I <y%w. ■ C « ■ E *("0*(d - "0 + ^ky ' ^(d) 

V 2 / m'=l \ 2 J ' 

n d+2 / r 



where Cq > is some constant. Now by Lemma 4.2, we obtain that 

n d+2 / r 



r(^)- ' V d-p! 
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Choose pi so small that C\C%p\ < \ and p± • 4a < \. Then take L so large that < \. 
This clearly implies 

Pi +2 

|^mi,m 2 l — y ( d+7 ) 

We now justify the induction hypothesis (29). Recall that our free parameters are , h^i 
and h^j . It is easy to check that if we set = h 2 } = = 0, then h mim2 = for any 
m l + m 2 > 2. Since L is fixed, and < \h^\ < 5, < \h^\ < 5, < \h$\ < 5 with 
sufficiently small 5, then the induction hypothesis is satisfied. A simple estimate on g gives 
that 

g{m) < (4a) m 
Thus the theorem is proven if one takes p = Aapi. 

As it is stated our solutions of (20) are determined by five parameters , h^i , , . 

However, it turns out that these parameters are not independent and o\ can be expressed 
through ( h$,h$,h$ ). Namely, let <3" (1 UMMV 2) (Y) be the solution of (20). Then 

G (*Uh<8,hU,hg>,vV) {Y) = G (l,*W(h<8 - 1) + l,*Wh$, *W(h§> - 1) + 1, aW) {Y) 

This equality is proved at the end of §6. We formulate now the final result concerning the 
existence of solutions of (21). 

Theorem 4.2. Let > 0, > and , ^4i\ ^4o be sufficiently small. Then there 
exists a solution of (20) which has the following form 

et^SMMVV,,^) = ex P {-^ (mp + mn} ■ 



a — ■ exp(-^-|r 3 | 2 ) J— • v^^'^^^v^y^y,). 

Here > h ^ ' h so ) is the solution of (21) with the given h!{ 2 \ , ■ 

As it was already mentioned the parameters a±, h^ 2 , , h^j , a 2 are not independent 
and actually the set of solutions depends on four independent parameters (see Lemma 6.2). 



From the estimate in Theorem 4.3 and from known asymptotic formulas for the Hermite 
ynomials it follows that the series givi 
Better estimates are also easily available. 



polynomials it follows that the series giving ij( ft -i2 ) '4i ) >4o ) ) converges for every Y = (Yi,Y 2 ) 
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5. The Linearization Near Fixed Point 



Denote = h 2 } = x^ 2 \ = x^ 3 \ Our fixed points have the following form 

r{a m, x m, x (v, x (v,*(V) _ ff(1) nvT J <r w (Y? + Yj) \ 
~ 2tt \ 2 J 

/ (30) 



Recall that #(- (1) ,* (1) ,* (2) ,* (3) ) = v /^#(W 1 W 2 \^)( v /^Wy 1) V^Y 2 ) and 
are described in §4. 

The strategy of the proof of the main result is based on the method of renormalization 

group. At the p-th step of our procedure, we consider an interval on the time axis = 

S (p \ S^ ] ] such that C S®. From our estimates it will follow that f| = [S-, S+] 

v 

is an interval of positive length. We want to find conditions under which g r (Y, s), s e 
have a representation 



(h^^^Y) + 5{ r \Y, s), H^^^^iY) + 5i r \Y, s), #>(y, 

where d[ r \ 8<f\ 5^ tend to zero as r — > oo. The renormalization is based on the crucial 
observation (see above) that for large p, the sum over p\ is a Riemannian integral sum for 
an integral over 7 changing from to 1. Let us write 

~, v s>- r +i, 1 f^ (1) (n 2 + ^ 2 ) a^Y, 2 \ f 2n\ f 2n\^ 
g r (Y,s)A r+1 (r x exp<^ y — — + ' ( ' 11 1 



2 2 J V^ (1) / W 2) / ' (31) 

= H^ 1) ' xW ' xi2) ' x ^{Y 1 ,Y 2 ) + 5 ir \ 1 ,Y,s) 

where 6^(7, Y,s) = |^(t, Y, s), 1 < j < 3} = ^(7, Y, s), 7 = T -. It is natural to 

consider the set of functions {#^(7, Y,s)} as a small perturbation of our fixed point (30). 
Recall that the third component of f/~(°" (1) . a;(1) . a;(2) > a;(3) ) j s zero because of incompressibility and 
5^ can be found from the incompressibility condition. Clearly, 

^)( 7 ,F, S ) = ^)(^ 7 ,F, S ) , 7<^r- 
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The formula for 5^ p+1 \l, Y, s) follows from (21): 



cxp 



(|yj 2 + |y 2 | 2 ) - 



"~'|2 , \xr\2\ ° \y\ 2 

2 1 31 

(2) 



<jW urn 



a 



(i) 



(i) 



2tt V 2tt 



r(2) 



K3 



2tt7 \/ 2tt7 2tt(1 - 7) V 2tt(1 - 7) 



exp 

a( 2 )|y 3 f 

2tt(1-7) 
+ 7 I(1 _ 7 ) 



27T7 



27T7 



i - 1 Yj Y[ rj (Y — Y' 
-(1 -7) 2 I — #1 



a/7 



#1 



y-y 



+ 



Y< 



a/7 



Ho 



2tt(1 -7) 



^2-^2 - /y-V 



#2 , 

\/7 V \/7 



y-y 

^7~ 



r(p+i) 



1-7, 



y 



+ 



(l _ 7) .(^( 7 .-:,) + ^ 



a/T=7 



f 7^(1 -7) (^=^(7,^^) + -7=^ +1) ( 



y 



d 6 Y', 



J = 1,2 



7 ' Vi ' 
y-y 



(32) 



We did not include in the last expression terms which are quadratic in 5 because in this 
section we consider only the linearized map. 

Another way to introduce the semi-group of linearized maps is the following. Take 9 > 
which later will tend to zero. Denote 7,- = (1 + 0)~\ j = 0, 1, 2, . . . Our semigroup will act 
on the space A of functions 5(7, Y) with values in C 3 such that 



1. for each 7, < 7 < 1, the function 5(7, Y) belongs to the Hilbert space L 2 = L 2 (i? 3 ) 

3 

of square-integrable functions with respect to the weight (j^Pj 2 exp{ — ^-j—}, Y = 
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2. As a function of 7 it is a continuous curve in this Hilbert space and max ||5(7, ^)||l 2 < 

0<7<1 

OO. 



Define the linearized map Lg corresponding to 6 as follows: 

1. for 7i+ i < 7 < 7j, j = 1,2,... 

L,(5(7,F)) = 5(7(1 + 0),y); 

2. for < 7 < 1 the function L 9 (5('-f, Y)) is given by the formula 

1 + 9 

L (6(>y,Y)) = 6 pl (l,Y,s) 

where p\ is found from the relation — = 7. 

V 

In other words at 7 = 1 we use (32) to find the new 5 ip+1 \l,Y,s). After that we apply 

1. 



It is easy to see that there exist the limits lim Lg = A f and the operators A 1 constitute 
a semi-group. For 7 < 1, t > such that 7e* < 1 

A t 5( 1 ,Y) = 5( 1 e t ,Y) 

Let ^4 be the infinitesimal generator of the semi-group A 1 . In §6 we study in more detail 
the spectrum and eigenfunctions of A. 

Lemma 5.1. The eigenfunctions of the group A 1 have the form 

<*( 7 ,y) =7 «$ a (y) 

where $ a is a function with values in C 3 satisfying (32). 

In more detail, if we take <5( 7 , Y) = 7 a $ a (F) and substitute it into the rhs of (32) we get 
in the Ihs 5 ip+1) (Y) = $ a (Y). 
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Proof. If 5(7, Y) is an eigenfunction then from the formula for A 1 



Let 7 -> 1. Then 



Lemma is proven. 



am(7, y) = 5(76-*, y) = e- Q M( 7 , y) 



5( e -*,y) = e- Qi $(y) = 7 °$(y) 



The space A is spanned by the eigenfunctions of {A 1 } in the sense that for any h G A 
we have the expansion 

H 1 ,Y)= c (a) 7 Q ^(y) 

aGspec.4 

The coefficients are found with the help of the eigenfunctions of the conjugate system 
{(A*) 1 }. The form of the conjugate semi-group and its eigenfunctions can be investigated 
using the described above discrete approximation. We do not dwell more on this. 

§6. The Spectrum of the Group of Linearized Maps 

In this section we show that the solutions of (21) studied in §4 have = 4 unsta- 
ble eigenvalues and = 6 neutral eigenvalues. Therefore in the renormalization group 
approach we consider 10- parameter families of initial conditions (see below). 

As was already mentioned, in the limit p — > 00 the linearized maps generate a semi- 
group of operators acting in the space A of functions /^(7, Y), < 7 < 1, Y e M 3 , j — 1,2 
which are continuous as functions of 7 in the Hilbert space L 2 . At 7 = 1, the functions 
/( J )(7,y) satisfy the boundary condition which follows from (32): 



d 7 







2tt7 V 2tt7 2vr(l - 7) V 2tt(1 - 7) 
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exp 



^(l^-y/p + l^-K,'! 2 ) aW\Y 3 -Ytf a^(\Yl\ 2 + \Y^\ 2 ) 



27T7 



2"7T7 



2tt(1 -7) 



q (2) |y 3 f 

2tt(1 - 7) 



a/7 



a/7 



a/7 



a/7 



y 



a/T^7 



+ 



1-7) 



f( i)/ y-y\ y 2 -y 2 ' 



a/7 



( 7 



a/7 



+ 



v/7 



/ (2) ( 7 



y-y 



+ 7 2 (I-7) 

Y' 

7r^ 



y 



( 7, 



y-y 



V 1 -7 ' V" a/7 

d 3 r, j = 1,2 



+ 



y 



a/T^7 



/( 2) 



7, 



a/7 

y-y 



(33) 



Denote by 1Z P the linear operator which transforms {5^ (7, y, s)} into {5( p+1 )(7, Here 
s is a parameter which plays no role in this section. As it was explained in §5, for each t 
there exists the limit lim IZf = A 1 so that the operators A 1 constitute a semi-group having 

p >oo 

an infinitesimal generator A = lim ^7-^ In our case .45(7, Y, s) = 7 dS ^' s ^ , < 7 < 1 and 
for 7 = 1 the function 5(1, Y, s) satisfies the boundary condition (33) in which 7(7, Y) = 

^ +1 )(i,y, s ). 

If a is an eigenvalue of A, then the corresponding eigenfunction has the form 7 Q( I > Q , )(J (i) i0 -(2) (Y) 
(see Lemma 5.1), where $ a CT (i) CT (2) (y) satisfies the equation (33) with 7(7, Y) = 7 a $ Q , (T (i) jCr (2) (Y). 
If dt(a) > (5ft(a) = 0) then the corresponding eigenvalue is called unstable (neutral). All 
other eigenvalues are called stable. The subspaces generated by unstable, neutral, stable 
eigenvalues are denoted by T^ u \ T^ n \ respectively. 

As before, for 3>^') (1) (2) (Y) the following scaling relation with respect to is valid: 



Therefore it is enough to consider the above equation (33) for o^ 1 ) = cr( 2 ) = 1. We again use 
the expansion over Hermite polynomials: 
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*S,iOn = *%>(Y) = J2 fi j) (m u m 2: m 3 )He mi (Y 1 )He m2 (Y 2 )He m3 (Y 3 ) 



mi, m,2, m3 



Here j takes values 1,2,3. Since in m 3 it is the usual convolution and H does not depend 
on y 3 , it is enough to look for solutions of (33) having the form / mi , m2 # m3 . Put /9 = a + ^y- 
and fp\mi,m,2) = fa\mi,m2)S m3 . We come to the linear system of recurrent relations 



fp\™u ™a) = £ J# m ,, +2/J K, m' 2 ) + (B 2 h^) (mi, m' 2 ) ) /^K, m' 2 ') 



+ 4^+ 2/3 • K, O ■ Wf) K, m' 2 ') 

+ 45U» +2/3 • fc (2) ("»i,m' 2 ) • (S 2 /W)K,mJ) 

+ ^ +3Am « " ((^tf'lK.mg + (<)KX))li«K,^ (34) 
+ • /JViX) ' (5i^ 0) )K,m' 2 ') 

+ J { Xw ■ • (s 2 ^))K',m' 2 o 

Introduce the vector 

4 d) = (if (0, d), 4 1} (1, d - 1) , . . . , ff ] (d, 0) 4 2) (0, d), (1, d - 1) , . . . , 4 2) (d, o)) T 

The vector /f contains all terms of degree d. The recurrent relation (34) can be written as 

where the vector contains terms of degree < d — 1. Also Cf € jg>2(<m)x2(c(+i) j g a ma ^ r j Xi 
Let Ci d \k,£) be its (Jfe, t) -entry. Then 
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cf(M) 



16d + 32/3 - 16 + 32k 



1 - 



(d + 2/3 + l)(d + 2/3 + 3)(d + 2(3 + 5) ' 

80d + 160/3 + 80 - 32A; 
(d + 2(3 + l)(d + 2/3 + 3)(d + 2/3 + 5) ' 

32(d + 2/3 - fc + 2) 
' (d + 2(3 + l)(d + 2/3 + 3)(d + 2/3 + 5)' 

32(fc — d — 2/3 — 1) 
(d + 2/3 + l)(d + 2/3 + 3)(d + 2/3 + 5)' 



I 0, 

Note that d + 2K(/3) > -I. 



i{l<k = ^<d + l 

tfd + 2<k = £<2d + 2 

H2<k<d + l,£ = d + k 

iid + 2<k<2d+l,e = k-d 
all other cases 



Lemma 6.1. Assume 9ft(/3) > 0. There exists an integer d* > 0, independent of (3, such that 
for all d > d* ; the matrix Cjf* is invertible. 

Proof. As d tends to infinity, Cjf* tends to the identity matrix if 9?(/3) > 0. A simple 
estimate on the diagonal and off-diagonal entries shows that for all (3 such that 9£(/3) > 
and sufficiently large d, the matrix cjf^ becomes diagonally dominant. This implies the 
existence of the needed d* and its independence of (3. 

A similar statement holds if we assume that 9?(/3) > —A for any given A < 0. We 
formulate it as the following lemma. 

Lemma 6.1'. For any A > 0, there exists an integer d*(A) > which depends only on A, 
such that for all d > d*(A) and all (3 with K(/3) > —A, the matrix cjf* is invertible. 

By Lemma 6.1, to find all eigen- values of A it amounts to solve the equation det(C^) = 0. 
The eigenvalue a is then found from the relation (3 = a + Let 

16 \ / f 32 



(d + 2/3 + 3)(d + 2/3 + 5)7 / \{d + 2/3 + l)(d + 2/3 + 3)(d + 2/3 + 5) 
Then a± is the eigen-value of the matrix E ^ 2 ( d + 1 ) x2 ( d + 1 ) given by: 
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f fc-1, 



C (d \k,£) = < 



itl<k = l<d+l 



2d + 2-k, ifd + 2<k = £<2d + 2 

d + 2-k, i£2<k<d+l,£ = d + k 

k-d-1 if d + 2<k< 2d + I, £ = k-d 
0, 



all other cases 

It is not difficult to find that the eigen- values of are and d+1 with algebraic multiplicity 
d + 2 and d respectively. Solve the equations a\ = or a\ = d + 1 and use the condition 
d + 23?(/3) > —1. The possible values of (3 are then given by 



3-d 
p = — - — or 



y/Vf-A-d 



d= 1, 2, 3, 



2 2 
This fact immediately gives the following lemma. 

Lemma 6.2. Let (C^)' 1 be th e inverse matrix of cf ] for d > d*{(3), where d*{(3) =3-2/5 
or y/Vf — 4 — 2(3 is an integer. Then there exists an absolute constant C 2 > such that for 
alld> d*(P) 



< Co. 



We now state our theorem about the properties of the solutions to the recurrent relation 
(34). 

Theorem 6.3. The only possible values of (3 for which (34) have nonzero solutions fp \m 1 ,m 2 ) 
is given by: 

The corresponding solutions /^'(mi,m 2 ) have the following property: 



a) (3 = (VT7 — 4 — m)/2. In this case /^'(m 1 ,m 2 ) = for any < m\ + m 2 < m. For 



d = m, we have 



ff\r,d-r) = -(d-r + l)ff\r-l,d-r + l) J r = l,2,---,d 



25 



f^\o,d), f^\d,0) are free parameters. /^(mi,m 2 ) for mi + m 2 > m+1 are uniquely 
determined if the values of the m + 2 free parameters f^\r, m — r),r = 0, 1, • • • , m, 
and f^\m,0) are specified. 



b) (3 — (3 — m)/2. In this case /^(mi,m 2 ) = for any < mi + m 2 < m. For d 



m 



we have f^\o,d) = f^\d,0) = 0, and 



(3 

4 1) (r,d-r)=4 2) (r-l,d-r + l), r = 1, • • • ,d 

are /ree parameters. /^(mi,m 2 ) /or mi + m 2 > m + 1 are uniquely determined if the 
values of the m free parameters f^\r, m — r), r = 1, • • • ,m are specified. 

In both case a) and b), the solutions /^'(mi,m 2 ) zs zero /or m x + m 2 = m + 1, m + 3, . . .. 
S'mce fp^ depends linearly on the free parameters, we have for some C3 > 0, < p < ^ 



fp J \m 1 ,m 2 ) 



4000 

pmi+m.2+2/3 

^ g 3 r ^ mi+m2+ W y , Vm!>0, m 2 > 0, J = 1,2. 



Proof. Property a) and b) are straightforward computations. From recurrent relation 
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(34), by parity it is obvious that /?(mi,m 2 ) = for mi + m 2 = m + 1. An easy induction 



shows that /^(mi,m 2 ) = for mi+m 2 = m + 3, m + 5, . . . We now prove the decay estimate. 
The strategy of the proof is the same as in theorem 4.3. From the proof of theorem 4.3, it 
is clear that by choosing the parameters ( ( 3 )) sufficiently small, we have 

\ h m\,m 2 \ < ^ fm ,^n+7\ Vm i > 0, m 2 > 0,mi +m 2 > 3, j = 1,2. 



Our induction hypothesis for /®(mi,m 2 ) is 

n m.i+m 2 +2/3 



Y ^ r»i+r» 2 +2/3+3 ^ 



< ~ / mi+m2+2/ 3 — Vm<mi+m 2 <d, j = 1,2. 



where d > L and d — m is an even number (note that /^(m 1 ,m 2 ) = for m x + m 2 = 
m + 1, m + 3, . . .). We assume that L is a sufficiently large number and will verify the 
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induction assumption for m < d < L later. Now for mi + m 2 = d, by lemma 6.2, we have 



d—m 



4 j) (m 1 ,m 2 ) <C 2 - ^(m' + l). 

m'=2 
d—m 

+C 2 • ^ (m' + 1) • 



m'=4 

+C 2 - I J, 



J (1) 

U m',m"+2P 



J (1) 

d m' ,m" +2(3 



2(m' + 1) 



P 



m'+3 



P 



m"+2/3 



p ^ m'+8 ^ y ^ "i"+2/3+3 ^ 



■ 2 



P 



m'+l 



P 



m"+2/3 



r(^) r( m " + f +3 ) 



a) 

2,d-2+2/3 



•4- 



d- 2+2/3 

r(^±f±i) 



d— m+1 



m'=3 
d— m— 1 

E 

m'=3 

' r (2) 

,d- 1+2/3 

d-3 

E 

m'=m— 1 
d— 3 

E 

m'=m+l 

' J (l) 
'd-1+2/3,1 

d-2 



+c 2 ■ K + !) 

m'=3 
d— m— 1 

+C 2 • ^ (m' + 1) 

m' 

+c 2 • I j{ 

d 

+C 2 - £ (m" + l). 

m'=m— 1 
d-3 

+C 2 - ]T (m" + l). 

m'= 

+c 2 -|j, 

+C 2 • ]T (m" + 1) 

m'=m 
d-A 

+c 2 ■ J2 K + 1) 

m'= 

+c 2 -U 



7(2) 

J m',m"+2P 



• 2 • 



P 



m'+2 



t(2) 

J m',m"+2P 



r(^ti) 

,m'+2 



(m" + 1) 



P 



m." +2/3+1 



p ^ to"+2/3+4 ^ 



P 



m"+2/3-l 



,d- 1+2/3 



p ^d+2/3+2-j 



J (1) 



7(1) 

J m'+2P,m" 



,d-2+2/3 



r(^+f±i) 



7(2) 

J m'+2P,m" 



2 • (m! + 1) 



P 



m'+2/3+l 



P 



,m"+2 



r(^) r(^ti) 



V+2/3-1 



y+2 



• 2 • 



p ( m< +2(3+3 ^ p ^ m"+7 ^ 



P 



m'+2/3 



P 



,m"+3 



p ^ m'+2/3+3 ^ p ^ m"+8 ^ 



(m" + 1) 



7(2) 

J m'+2f3,m' 



■ 2 • 



P 



m'+2/3 



m—m 
' 7(2) 



r(=i^±3) r(^te) 



'd-2+2/3,2 



d-2+2/3 



r(«+f±i) 
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P ^ ( d ^ 8 P 3 ^ 8 p 8 

-°2 r /d+2|+3N lZ^ rf + 2/ 3 + 5 + + 5 p 2 -(rf + 2/? + 5) 

V 2 / \m =2 rrt '= 4 

"y^ 1 8 P 3 V + 1) 'y^ 8 P 8 

m'=3 m'=3 



r 

d— 2 , „ o d-A 



^ 8p 3 2(m' + l) ^4 8p 16 

+ ^ d + 2/3 + 5'd+2/3 + 3 + + ^ d + 2/3 + 5 + p 2 (rf + 2/3 + 5) 

m'=m-l ^ ^ m'=m+l ^ r v r > 



v 16p 3 v 8p 16 



^ d + 2(3 + 5 ^ d + 2(3 + 5 p 2 (d + 2(3 + 5) 



m'=m 



d+2f3 d+2/3 

<C 2 - "^ -2000p < ' 



where in the second last inequality above we have used the fact that d> L and L is sufficiently 
large such that d/(d + 2(3 + 3) < 2. It is clear that it suffices for us to take L = 2m. To check 
the inductive assumption for m < d < 2m, we recall that f^\mi,m 2 ) depends linearly on 
several free parameters. If we let them be sufficiently small, then it is clear that the inductive 
assumption is satisfied for m < d < 2m. Our theorem is proved. 

We now formulate our main theorem about the spectrum of the linearized operator. 
Theorem 6.4. The spectrum of the operator A consists of the following eigen-values 

spec(^l) = |l,i 0,AS\Ag>,m>l} . 
where = -f , = y^=m , m > i. 

The first eigen-values have multiplicities V\ — 1, v\ — 3, u = 6. The eigen-values Xm\ 
Am correspond to the stable part of the spectrum and also have finite multiplicities given by: 



_ (m+3)(m+4) _ m(m+5) 
~~ 2 ' U \ { - 2) ~ 2 ' 



For each a G spec (^4) , the eigenf unctions fa\m\, m 2 , m 3 ) have the following property: 

a) fa\mi,m 2 ,m 3 ) is compactly supported in the m 3 variable, i.e., there exists an integer 
m l = m a( a ) su °h that 

f^ J \m 1 ,m 2 , m 3 ) = if m 3 > m* 3 
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b) /a '(mi,m2,m 3 ) decays faster than exponentially, more precisely, there exist constants 
C3 = C^{a) > and < p < such that 

|/^' ) (m 1 ,m 2 ,m 3 )| < ^ mi+ma+m3+2a+3 ^ Vm 1 ,m 2 ,m 3 >0 

T/ie system of eigenfunctions is complete in the following sense. Let r (s) 6e the stable 
linear subspace of A generated by all eigenfunctions with 3?(A) < ; be the unstable sub- 
space generated by all eigenfunctions with eigenvalues A > ; and be the neutral subspace 
generated by all eigenfunctions with eigenvalue A = 0. Then dimT^ = 4, dimT^ = 6 and 



a = + + . 

Proof. By Lemma 6.1, we only need to examine P for which det(C^ ) ) = 0. From 
previous arguments, we have that for d > 1, (3 — — ^ or ^~ A ~ d , We discuss the spectrum 
separately in the following three cases. 

1° unstable spectrum: a = 1,1/2. 

a) a — 1. Since f3 = a + the only possibility is that (3 = 1, d = 1 and m 3 = 0. 
The eigenspace is one- dimensional with f$ = = / ( }J = = 0, /[JJ = / ( (2) 



010 



is a free parameter and the remaining part of all higher degree terms ( m2 with 
mi + m 2 > 2 ) is uniquely determined once we specify / 



(i) 
100 • 



b) a = 1/2. Possible cases are m 3 = 0, (3 = 1/2, d — 0, 2 or m 3 = 1, /3 = 1, d — 1. In the 



first case we have /^^^o = for mi + m 2 < 1, = fj®, f$ = /gj are two free 

fCj) 



parameters, all other terms of higher degee ( f^ i m2 with m x +m 2 > 3 ) are uniquely 
determined once we specify the above four parameters. In the second case we have 
/001 = /001 = /on = /101 = 0> /101 = /on is a free parameter and the remaining part of 
all higher degree terms ( f^ 1Tm i with m x + m 2 > 2 ) is uniquely determined once we 
specify f[l\. Putting two cases together, we see that the dimension of the eigenspace 
is 3. 
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This gives dimr (M) = 4. 

2° neutral spectrum: Here we have a — 0, and three cases. 

a) m 3 = 2. Then (3=1. The eigenspace is one-dimensional with /qq 2 = /002 = /012 = 
/102 = 0' /102 = /012 i s a f ree parameter and the remaining part of all higher degree 
terms ( m2 2 with mi + m 2 > 2 ) is uniquely determined once we specify f$ 2 - This 
eigenvector is connected with -^^j which corresponds to the variation of the parameter 
(j( 2 ) of the fixed point. 

b) m 3 = 1. Then (3 = 1/2. We have /£> >IBa>1 = for m x + m 2 < 1, /<}} = / ( 2 2 j, /£} = /£> 
are two free parameters, all other terms of higher degee ( /^ m2]1 with mi +m 2 > 3) are 
uniquely determined once we specify the above two parameters. Clearly the eigenspace 
is two-dimensional. This eigenspace does not correspond to any change of parameters 
of the fixed point. 

c) m 3 = 0. Then (3 = 0. We have f^ m0 = for m x + m 2 < 2, /£ = /£ = 0, 
fill = io3Cb /210 = A ( 2o> /So = /210 are three free parameters. All other terms of higher 
degee ( m2 with mi + m 2 > 4 ) are uniquely determined once we specify the above 
three parameters. This eigenspace corresponds to (^ry, afW)- 

Putting all three cases together, we see that dimT^ = 6. 
3° stable spectrum: $l(ct) < 0. 
There are two cases. 



Case 1 : a = — y, m > 1. Recall that (3 = a + and m 3 satisfies < m 3 < m + 2. By 
theorem 6.3, for each such (3, the number of free parameters is 3 — 2(3. Then the 
total multiplicity u a is given by 

m 3 =0 
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Case 2 : a = ^~ 4 ~ m , m > 1. /3 = a + and m 3 satisfies < m 3 < m — 1. By theorem 
6.3, we have 

E, ^ m(m + 5) 
(m - m 3 + 2) = 

rri3=0 

It follows easily that the eigenfunctions fa\m 1 ,m2,m 3 ) is compactly supported in the m 3 
variable. By theorem 6.3, the decay estimate on fa\mi,m 2 ,m 3 ) is obvious. 



It turns out that the eigenvector corresponding to -q^jj is in the eigenspace spanned by 
the eigenvectors (gjry, g^2j, &fW)- More precisely we have the following: 

Lemma 6.3. Let t x = x w - I, t 2 = x( 2) , t 3 = x (3) - 1. Then 



where q(° {1) > x(1) > x(2) > x(i) >° (2) ) i s defined in (30). The function G satisfies the following scaling 



relation: 



t( 2 ) 



(35) 



Proof. Let /^'^ 2)0 correspond to the eigenvector ^tj, then a simple calculation shows 
that 



(?) 



,0") 

mi,m2-2 ' 



^ /mi,m2,0' fmi,m2,0 and fmi,m 2 fl correspond to the eigenvectors and 
respectively, then clearly we have 

f Ci).0 _ ( (1) _ f (J),l , (2) f (J).2 , ( (3) _ i \ f (j).3 
Jmi,m2,0 ~~ ^ x V Jmi,m2,0 ~ x Jmi,m2,0 ~ ^ x V Jmi,ro 2 ,0 

This immediately gives 



r(1) ^-(- (1) - 1 )^T) - (2) CT-(* (3) -l) ' 
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Regarding this as a transport equation in the variables (a^\ti, t 2 ,t 3 ), we can easily find that 
G satisfies the scaling (35). Lemma is proved. 

This lemma actually shows in what sense the parameters x^ 2 \ x^ are dependent. 

As was shown in §4, we have the five-parameter family of fixed points q( u{1) < x{1) > x{2) < xii) ^ 2) ) . 
We use the notation n = (a^\ x^\ x^ 2 \ x^\ a^) and write instead of (^( cr(1) . a;(1) . a;(2) . a;( ' :!) . o " (2) ) 
The spectrum of the linearization of the equation for the fixed point does not depend on it 
(see §5) and has = 4 unstable eigenvectors $J u) (Yi, Y 2 , Y 3 ), l<j< £ {u) = 4 and £ {n) = 6 
neutral eigenvectors $$°(Yi, Y 2 , Y 3 ), 1 < f < £^ = 6. 

§7. The Choice of Initial Conditions and the Initial Part of the Inductive 

Procedure 

The equation (21) for the fixed point which was derived in §3 is non- typical from the 
point of view of the renormalization group theory because it contains the integration over 7, 
< 7 < 1. On the other hand, since we consider the Cauchy problem for (1) we are given 
only the initial condition v (k, 0) which produces through the recurrent relations (4), (5), (6) 

or (4'), (5'), (6') the whole set of functions h r (k, s) or g r (k, s). For large p and r < p they 

r 

can be considered as depending on 7 = - and our procedure is organized in such a way that 

P 

for 7 which are away from zero g r are close to their limits. Therefore the initial part of our 
process should be discussed in more detail. This is done in this section. 

We take k^ which will be assumed to be sufficiently large, introduce the neighborhood 

A 1 = {k : \k-K^ \ < D^kWlnkW} 

where = (0,0, k^) and D x is also sufficiently large. Our initial conditions will be zero 
outside Ai. Inside A x they have the form 




6 

bf^f^ Y 2 , Y 3 ) + $(Y"i, Y 2 , Y 3 ; b^, &<»>) 

j'=i 
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In this expression k = fc (0) + Vk^Y, H^(Y U Y 2 ) = (h[°\Yi, Y 2 ), H { 2 ] '(Y, Y 2 ), 0) is the fixed 
point of our renormalization group (see §4) corresponding to the parameters = = 
1, xi — x 2 — x 3 — 0. Also $^ are unstable and neutral eigen-functions for 

described in §6, and are our main parameters, —p\ < b^\ < p x where p\ is 
another constant which depends on A;( ). Its value will also be specified later. Each function 
$(yi,y 2 ,*3;& (u) ,& (n) ), = {bf ] }, = is small in the sense that they satisfy 

sup|$(Yi,Y 2 ,Y 3 ;&M <D 2 , 

Y,b 

sup ||$(Y 1? Y 2 ,Y 3 ; 6^,6^) - $(Yi, Y 2 , Y 3 ; 6 (n) )|| < £) 2 (|6^ -6^| + |& (n) -& (n) |). 

Due to the presence of 6^ M \ 6^ n \ we have / = + = 10-parameter families of initial 
conditions, due to the presence of <3> we have an open set in the space of such families. 

Let 

A r = {k: \k- r/t (0) | < D^rk^lnrkW} 

and the variable Y be such that k = rn^ + V rk^Y. Assume that for r < p, \Y\ < 
DiVlnrfcW 

/i r (r/t (0) + Vri^Y, s) = Z p (s)k r -\s)r~g r {Y, s) 

and 



g r (Y,s) =-exp | _ j — exp |_ 



#} 0) (Yi, Y 2 ) + ^(Y, Y 2 , Y 3 ), ^ 0) (Y, Y 2 ) + 4 r) (Y, Y 2 , Y 3 ), 



where in view of incompressibility 

H[ 0) Y 1 + Hf ] Y 2 + F (r) = (36) 

We shall derive a system of recurrent relations for Z p (s) and A p (s) for p < p . All 5 ? - r * ) 
will be considered as remainder terms. 

Outside A r we assume that 
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| M ™(o) + v ^y, s) |<_^ 
where Ai is another constant which depends on C\. 

Returning back to (6) take the term with some pi, p 2 , P1+P2 — P and introduce the new 
variable of integration Y' where k' = p 2 n^ + y/pk^Y'. Introduce also the variables 9 ly 9 2 , 

0<9 1 < ( Plk My, 0<9 2 < (p 2 k^Y where 8l = s (l - , s 2 = s (l - . 
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Then from (6) 



h p (p^ + y/pk^Y, s) = Z p+1 (s)A p p+1 (s)p~g p (Y, s) 

(( P -i)fc(°)) 2 

= (pk^)ii [ d6 2 [ exp [-9 2 \kW + - Y ' | 2 







&((y - y)v£fcw, o), «<°.°> + -^=) p K(0 , 0)+ ^_g P -i (V\/A 8(i - -4— )) d*r+ 

Vp^/ Vpfcw V V p- 1 (( P ~i)fc(°>) 2 / 

(pifc(°>) 2 (p 2 fc(°») 2 3 

1 (pkW)*pip 2 f f f ( 1 V f (n-n^+^-v^+re.-^') 2 







p(pikwy(p 2 km 2 j L j J \2icj *\ 2 7 

^7 R3 



Z p (s(l — )) • A^ 1 ~ 1 (s(l — )) • ZJs{l — )) • K^\s{l — ))• 



1 y_\ «)W^U| gl | K (Q,Q) + Z^Lpi exp j- g2 [^°) + r ' I 2 I + 



— exp < > exp < —pi k H ; /"exp. 

27T / [ 2(1-7) J I -r^/pk(°) J ^ (l-7)\/pfeW 



((p-l)fc(O)) 2 



y " y |2 



R3 

(37) 

Here 7 = — and /t(°>°) = (0,0, 1). Now we shall modify (1371) for pi > 1, p 2 > 1 similar 

p 

to what we did in §3. Later we discuss the terms with p 1 = 1 and p 2 = 1. The modification 
consists of four steps. 

Step 1. All terms s(l- (pi f(o ))2 ) ,s (l - (p J(o ))2 ) are replaced by s. 
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Step 2. Write 



(pfc(°))5 



(pifc(°)) 2 (p2fc (0) ) 2 (fc(°)) 2 7(l - 7) 
Step 3. Consider the inner product 



Up to remainders and from (J36|) it equals to 



3 

1 N 5 
2^ 



exp 



(Y 1 - Yjf + (y 2 - yp 2 + (y 3 - K 

27 



'A2 



(0) ^1 - *i *2 ~ ^ V , H-(O) f Yl ~ F l ~~ F 2 \ y + J__p(pi) 



v^7 



>/7 v 7 ^ 



y-y 



„ 3 
1 \ 5 

2^ 



exp 



(y x - y/) 2 + (y 2 - y 2 ') 2 + (y 3 - r 3 ') 2 



(0) (Y 1 -Y{ Y 2 -Y 2 ' 



v/7 



7 



2 7 

(o) f Yi — Y{ Y 2 — Y 2 ' 



Y, 



(0) fy - y{ y 2 - y 2 '^ y 1 - y; ^ (0) fy - y[ y 2 - y 2 ^ y 2 - y 



7 



7 



3 

1 \ 5 
2^ 



exp 



(y - y/) 2 + (y 2 - yp 2 + (y 3 - y 
27 



'A2 



7-1 



7 



(0) f jfr - Y[ Y 2 -Y> \ Y 1 - Yj ^ h(0) (Y x - Yj Y 2 -Y> \ Y 2 -Y> 



x/7 



7 



V7 



Hi 



(0) / ^1 — Y[ Y 2 — Y 2 ' \ Y 1 



Vi Vi J V 1 -7 



H. 



\fi V7 y V7 

(o) / — ^1 ^2 — ^2 \ ^2 



+ 



V 7 ^ y/l ) V 1 -7. 



Let us stress again that Hj(Y, s) depend only on Yi, Y 2 and s. With respect to Y3 we have 
the usual convolution. 



Step 4. Replace the projection operator by the identity operator. It is not the reduction 
to the Burgers system because the incompressibility condition is preserved. 
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Now we shall modify the first and the last terms in (|37|). For the first one we can write 

((p-l)fc(0))2 



R3 

«(o) + (y _ y') V^W| 2 ) / v („(o) + (y _ y') vSfeW, 0), K (°' > + -^=\ • (38) 

•P (ou y fif p _! (Y\r^,s(l )V r 



• exp < —s 



The factor (p — 1) comes from the inductive assumption concerning h p -i. As before, we 
replace exp < — ^ 2 |^ ' ^ H — , y ' | 2 > by exp{— 6 2 }, P m , y by the identity operator 

and (^J^,s(l - ((p-fffcwja )) b y 9 P -i( Y '\J^i, s). All corrections are included in 
the remainder terms. 

For the Gaussian term in t> (k^ ^ + (Y — Y') y/pk(°\ 0) we can write $ exp j ^ y ~ y I p | . 

(27r)2 " L J 

This shows that typically Y — Y' = 0(A=). For the third component of v(k^ + (Y — 
using the incompressibility condition we can write 

(«(o) + (y _ y') , 0) = 



p\Y-Y 



n2 



■ exp 

For the inner product in (}3*Hj) 

y \ f p|r-r'| 2 

2 



H — -j= ) = exp 

\ ypfc(°) / 

• "#< 0) ((y - y')v$m + #(°)((y _ Y')^)Y 2 - 
- VP (<Xi - Y;)VPH?\(Y - Y')y/p) + (Y 2 - Y 2 r )^fp~Hf\(Y - Y')^p)) + 0(-±= 

The expression in the square brackets grows as y/p and therefore 
y/pk® ( v ( k (°) + (Y- Y') V^W, 0), K m + -L= 

\ v/pfeC ) 
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can be replaced by 



-y/p (Y 1 - Y{)^H?\{Y - Y')^p) + (Y 2 - Yi)y/pH?"((Y - Y')^p) 



(0), 



Further, 



exp{-s|«; (0) + (y - Y')^pW)\ 2 } = exp{-s|fc(°)| 2 }• 



•exp{-2sA; (0) (K (0 ' 0) , (Y - Y')^Vk^)} exp{-s\Y - Y'\ 2 pk^} 

The first factor takes values 0(1), the others can be written as 1 + 0(-t=). The main order 
of magnitude of (jHBj) takes the form 



pexp{-s(k^) 2 } 



(p-l) l 



(y - y/)^{ 0) ((^ - r')v$) + - nOv^HCr - r )v^) 



(0), 



+ 



H[°\(Y - Y')^p)Y x + H?>((Y - Y')^)Y 2 



, 3 

2^ 6XP 



\Y-Y'\ 2 p 



3 

1 N 5 
2^ 



exp 



2(P-1) J \ Vp-1 



A similar expression can be written for the last term in (}3Tj) . Remark that due to our choice 
of the interval S w the product s(k {0) ) 2 = 0(1). 



Now we derive the recurrent formula for Z p (s) and A p (s). Since our special fixed point 
is a Hermite polynomial of first degree, the convolution of over Y' (see (JH7J)) gives 

us simply the product of and the Gaussian term and a polynomial in 7. The function 
and the Gaussian term can then be taken out of the summation in 7 and this gives us 

the following recurrent system for Z p (s) and A p (s): 

Z p+1 (s)k p p+l (s) 

= E ~ v ITkyz ■ ^ - + 4 ) • ■ A » ■ C 1 - e- ( - fc(0))2 ) • (1 - e-toW) 

Pl+P2=P 

(39) 
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where the factor (67 2 — IO7 + 4) comes from the convolution of with itself. Now if we 
take Z p (s) = -i(k^ f and write = 1 + then we have 



-s(p 2 fc (0) ) 2 



) 



then it is not difficult to see that there exists bounded (with an bound independent of 
p) such that the equality holds. It is an elementary fact that the limit 



AM 



itaA P+l(s) =A 1 n(i + %i)' 

k=l v 7 



exists. 

Now we discuss the behavior of all remainders for p < [k^) X2 . 

By we denote the eigen- vectors of the linearized renormalization group corre- 

sponding to the fixed point H^°\ For each p we make the following inductive assumption for 

5^(Y,s), r <p: 

s"\y, .,) = ± (6« + /?<*>) *f + 1 (»# + + -r = r^r 

where 6^ = (p — 1)^6^7^' } , 6^ = 6^, is a function which belongs to the stable 
subspace of the linearized renormalization group, 7 = 

As we go from p — 1 to p, the variable 7 = ^-j- is replaced by 7' = ^ = 7 ■ Therefore 



K P — 1 

In the same way for the neutral eigen-functions we have 

because a,/ = 0. In the same way one can transform $( si ).The coefficients Z^]. are small 

•* Ji' J v 

compared to the first term . An important conclusion is that the projections to the unstable 
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directions increase, projections to the neutral directions remain the same and projections 
to the stable directions decrease. As was already said, in the case of unstable and neutral 
directions the term containing or is the main term. 

Now we discuss the form of S^(Y : s). It is the sum of three types of terms. 



ai). The term which depends linearly on all S^(Y : s). Especially important is the part 
which contains we were to have and be in the limiting 

regime then the integral will give p aj bj^ ^1 + $^ — (p+ l^^b^^f 1 since 
7' = 1. However, are slightly different from H^. Therefor we shall have a small 
correction which is included in all /3j£, f3^ p and in <& p st \ We denote it as (3^\, /3^\, 
®p*i- The we have terms which are linear functions of all /3j£, flftl and $l st \ They 
will give us $l,P%l*<$. 

a 2 ). The term which is the sum of all quadratic expressions depending on 5^ P2 \ We 
expand it using our basis of and all stable eigen- vectors. 

a 3 ). The term which contains all corrections which arise during the four steps described 
above. We also expand it in the same way as in 02)). 



The sum of all terms gives p^J , (3^,, 

We use this procedure till p = p = (k^) X2 . The procedure for p > p is discussed in §9. 



8. The List of Remainders and Their Estimates 



In the beginning of §7 we described 10-parameter families of initial conditions which we 
consider in this paper. We mentioned above that for each p we have an interval = 



S+* on the time axis. Actually these intervals are changing when p = p n = (1 + e) n 



where e > is a constant. Therefore we shall write 
will be no confusion. 



and hope that there 



In this and the next section we consider p > (k^) X2 . Each function g r (Y, s), 3 < r < p, 
has the following representation: 
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in the domain \Y\ < CiVlnrfcW, Y = (Y u Y 2 , Y 3 ) e TZ 3 

a™ 



g r (Y,s) = A r 1 • r • exp 



+ |i' 2 | 2 )} ■ \/^ ex p{-^W} ■ (/f""(y„n) + *"•'(!',»)); 



in the domain \Y\ > C\^hi{rk^) : 



2tt " [ 2 Vl " 1 " 7 J V 27r 

• ex P |-^|y 3 | 2 } • ^(y,^) + 5 w (y, s ) 



A3-1 



for some constant A3 > 0. We use the formula (7) to get g^(Y : s). New remainders appear 
in one of the following ways. 



Type I. The recurrent relation (7) does not coincide with the equation for the fixed point 
and actually is some perturbation of this equation. The difference produces some 
remainders which tend to zero as p — > 00. 

Type II. For the limiting equation all eigen-vectors in the linear approximation are multi- 
plied by some constant. In the equation (7) it is no longer true and the difference 
generate some remainders, (see also §9). 

Type III. The remainders which are quadratic functions of all previous remainders. 
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§8A. The Remainders of Type I. 

We call the domain A the set {\Y\ < Di^\n{rk^)} and the domain B the set 
{|y| > Di^/ln(rk^°y)}. The estimates will be done separately in each domain. We in- 
clude the first, the second and the last two terms in (7) in the remainders. We shall estimate 
only the first one, the others are estimated in the same way. 

Domain A : We have 

$\y,s) = (p+ i)i . -L . JdhJ <v + ^f- v^TT, o) ; , 

M.3 



d 3 Y' 



Here b means the collection of all parameters of v(k; 0). The the main contribution to the 
integral comes from Y — Y' = O ^ v /^+t ) • ^ n ^ ms domain in the main order of magnitude 

(^(o) + IjllL V ^TT, 0; b), yfikW) = 0(1) 

V s 

Assuming that v(k^ + \Jp + 1, 0; 6) is different iable w.r.t the first three variables we 
see that the inner product 

{v{k (o) + ^lT^Xo; a) >v ^fc(°) + -^=) 
V s vT + 1 

is of order 0(1). For g p we can write using our inductive assumptions 



9 P 



^ (2) (|y 3 | 



<^pi -.--.> • w" (nh 4) 
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Also 



y/sk^p + Y'y/p+1 



exp 



P 



and in the main order of magnitude the integration over 62 does not depend on Y'. Thus we 
can write 



\{3W(Y,s)\ < -p- exp^-^dm 2 + |F 2 | 2i 



(i) 



exp 



a 



(2) 



\Y*\ 



Da 
p 



(40) 



Here and later various constants whose exact values play no role in the arguments will be 
denoted by the letter D with indices. Since in the expression for g p+ i we have the factors 
A p ■ (p + 1) ■ expj-^dril 2 + |F 2 | 2 } ■ s£ • exp{-^|r 3 | 2 }, the estimate 



()40|) shows that \f3p (Y,s)\ is relatively smaller than g p+ i with an order O(^). This is good 
enough for our purposes. We did not discuss the errors which follow from the fact that the 
expressions in the previous formulas depend on 82- 



Domain B : The smallness of (3^ ( Y, s) in this case follows easily from several inequalities 
and arguments. 



1°: 
2°: 
3°: 



\Y\ < D 4 y/pkW because \k\ < D 5 pk {0 \ 



\Y — Y'\ < DqV because v(k, 0; b) has a compact support. 
2s. 



If \Y — Y'\ < 



If \Y - Y'\ > 



2s + 

VP 



then 



exp 



then 



fe(o) + (y - y') VpTT 



< 1 



exp 



< exp < — - 



Y-Y 



f} 
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4°: If \Y'\ > D 7 ^/p then 



cxp . 



yfik^p + r , v / pTT > < exp{-C 8 # 2 } 



5°: If \Y'\ < D 7 ^/p then 



fc (0) p + Y' y/p+1 



< 1 



6°: We have 



exp 



(i) 



r(2) 



mv + inn - v \y £ 



/|2 



= exp {- ^ (|y x - - YT)| 2 + |r 2 - (y 2 - r 2 ')| 2 ) 



r(2) 



in-^s-n' 



/\|2 



a 



(i) 



(7 



(2) 



= ex P - — (|Fi| 2 + |r 2 | 2 ) - — (|y 3 | 



• exp {a« {Y X {Y X - Y[) + Y 2 (Y 2 - F 2 ')) + (7® F 3 (>3 - K,') 



( 2 ) ^ 

(|r 1 -r 1 '| 2 + |y 2 -r 2 '| 2 ) - ^- 1^ 3 - v^| 2 | 



2s 



If |y-y'| < — then 



exp {,w (y^ - r;) + r 2 (y 2 - r 2 ')) + ^ 2 )( s ) y 3 (y 3 - KQ 



r(l) 



(l^-y/l 2 + \y 2 -y>\ 2 ) 



r(2) 



\Y 3 -Ytf \ < C 8 . 
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2s+ 

If \Y — Y'\ > then we have an integral of the function which is the product of some 

Gaussian factor and \H.^(Y)\. Direct estimate shows as before that in this case 



3 

p2 



which is also good for us. 

In the same way one can estimate terms with relatively small p\ and p — pi (i.e., p± < y/p 
or pi > p — y/p. The remainders will be of order — 1= • K The next set of remainders comes 
from splitting the integration over and Y' (see (7) and beginning of §3). We may assume 
that pi > y/p or pi < p — y/p because other terms were estimated before. Put 

p? v% 

|p+i(y,s) = t{ P + 1)5 V / de 1 / do 2 ■ — 2 



P1+P2 = P+1 
P1.P2 > VP 







/ 



• exp (-^ I vW> + .L-T | 2 - 2 | ViifcW + -^H*! 
I VP + 1 • 7 VP+I7 J 

Using the inductive assumption we can rewrite the last expression as follows: 



~g p+1 (Y,s) = i(p + l) J2 \ d ^\ 



Pi P2 

d9o 



P1+P2=P+! q 
P1,P2>Vp 



A-* ■ A P2_1 ■ _L_ ■ ^exp (-^ 

7(1-7) P+l I 2 7 
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2 7 2 (l- 7 ) 

-^•JS}.^<««(-v..(i-s)). 

As was explained before, due to incompressibility in the Domain A , the inner product 

(y-y.;.(i-|)),^ + ^> 

takes values 0(-j=) because the first two components of the vector y/sk^ + are of order 
0(-j=). Therefore the product 

Vp(h^ (V-y', s (i-|)) ,yTskW + ^=) 

takes values of order 0(1). 

The remainder can be written in the following form: 



p'i v\ 



P1+P2 =P+1 
P1.P2 > VP 







A pi-i . A P 2 -i / exp <^ - 



2 7 ' 2 7 2(1 -7) 2(1 



5} <««(,_ J)). 
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exp j-^*^ + ^^| 2 - 6? 2 |v^fcW + r 12 



V^(l-7) 



P1+P2 = J>+ 1 
P1.P2 > 1 



/exp[ 



7 (1 



p? pi 

1 1 f f 

- I exp {—9 1 s} dOi I exp{— 9 2 s} d6 2 

-V P J J 



a W (|y _ y/p + |y 2 _ y 2 f ) aW(\Y 3 - Ytf) 



27 



27 



^(i) (jy/p + y 2 f) |y 3 'p 



2(1-7) 2(1-7) 



• P 2 . (^(pi) (y _ y'), ^/ifcCo) + — ) P y W&»> (Y', s) d 3 Y' . 

VP y/sk {0) H 

We did not include the factor A p_1 • p because it is a part of the inductive assumption. This 
remainder is estimated in the following way. 

First we consider 



Ri = 



v^ (0) + 



Y-Y' 



s + 



v^ (0) + 



Y' 



As before, consider the domain where 

\Y-Y'\ < D 9 sj\n (pjfc(o)), |y'| < Dio ^/ln (pfc(°)) . 



We write 
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Ri 

In the Domain A 



\Y-Y'\ 2 \Y'\ 2 „ (\Y-Y'\ Y'\ 
V ■ 7i P • 72 



V^7l >/p( 1- 7) 



l^il < 



C 12 m(p£;( )) 



Therefore 



i? 2 = exp < -Q x 



exp{— #is} • exp{— 



v^7 2 



= ex V {-(9 1 + 9 2 )s)} 



exp < — #1 



Y — Y' 



exp < — 2 



V** (0) + 



v^a-7) 



and in the Domain A 



. 3 1 ■ Ci 3 #2 Inp 
I -R2 1 < exp{-(0i + 6*2)5} I h 



\/P7 v / P( 1 ~ 7) 



This shows that in the Domain A we can replace the exponent 



exp j-^lv^ + ^~\ 2 ~ 9 2 \yfthf® + Y ' 12 



by exp{ — (6 ) i+6 ) 2)s(/c^) 2 } and the remainder will be not more than^|p. This is enough for 
our purposes. In the Domain B the estimates are similar because again the main contribution 
to the integral comes from \Y — Y'\ < Dgy/ln p, \Y'\ < Dio\/\n p. In other words, in the 
Domain B we can replace the product of the Gaussian factors and by 



exp 



{-l* {1) (\Yi-Y;\*+\Y 2 -Y 2 r) - ± a W\Y 3 -n\ 
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\^H\Y{\ 2 + \Ytf)-\° {2 \\Yi\ 2 ))- 



This is also enough for our purpose. 

The next remainder of Type I comes from the difference between the sum over 7 and the 
corresponding integral. The remainder (3p (Y,s) is the difference between the sum 



v / 7V / Ft) • - • / expi- 



a W(\ Yl -Y{\ 2 + \Y 2 -Ytf) 



P1+P2 = 
P1.P2 > VP 



2 7 

(7 (2) ( |y 3 _y 3 /| ) 2 ^Dfly^ + |y 2 ']2) a (2)|y 3 /| 



2 7 2(1 -7) 2(1 -7) 



- 

1 v "' / ! J 2 . (w(pi)((y-y')),ViA; (0) + 5 



2tt7/ \2 7r (l ~ 7)/ ' 



and the corresponding integral over 7 from to 1. It is easy to check that this difference is 
not more than 

§8B. The Remainders of Type II and III 



All remainders of Type II appear because we use the sums (over pi) instead of the 
integrals . The functions H { ^^' ) are defined for all 7. We use a linear interpolation to 
define 5(7, Y, s) for all 7. From our inductive assumptions it follows that \8 p (i, Y,s)\ < 
Therefore, the remainders which follow from the difference between the sum and the integral 
also satisfy this estimate. 

It remains to consider quadratic expressions of 5 P (7, Y, s). The Gaussian density is present 
in all these expressions. Therefore, all the remainders are not more than 
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9. Final Steps in the Proof of the Main Result 



In this section we consider our procedure for p > p . Introduce the sequence p n , p n = 
(1 + e)p n -i = (1 + e) n p , where e > is small (see below). These are the values of p when 
we make the renormalization of our parameters. For p ^ p n , no renormalization is done. 

In §7 we explained the choice of our fixed point H (°\ The corresponding eigen- functions 
are denoted by $^ and Also we have eigen-functions of the stable part of the spectrum. 
Consider p, p m < p < p m+ \. By induction we assume that we have an interval on the time 



axis 



and s G 



, r < p, so that 



g r (Y, s) =A r ~ 1 ■ r ■ (H {0 \Y) + 5 {r) (Y, s))- 



^r exp 1 2 — rV^r exp 

If 7 = —^—r then 

' p— l 

5 {r \Y, s) = f: (&$ + ffi) ^^(Y) + £ (#i + * { r\Y) + *W{Y, 7 ). 

here (3^},f3^\ are small corrections to the main terms b^, , $r S ^ can be written as a 
series w.r.t. the stable eigen-functions. (see Appendix II). 

At one step of our procedure p — 1 is replaced by p, 7 is replaced by 7' = 7 • ^— ^ and ^ ) 

(u) («) 

is replaced by (l + ^f 3 ■ (7'P , b$ + 0$ is replaced by (5<$ + /#>) (l + . 
During the whole interval p m < p < p m +i the variable b^ m acquires the factor 



n 1+ 



(«) 

p — i j 

Pm<P<Pm + l X 7 

A similar statement holds for the stable part of the spectrum. The neutral part remains the 
same since = 0. 

Now we shall discuss S^(Y, s) using (7). As in §7 S^(Y, s) consists of three parts. 

Part I. In all 5^ the main term is the one which contains our basic parameters b^\ bp\ We 
consider terms in (7) which are linear in b^\ b^\ As it follows from the definition 
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of the linearized group and its spectrum we get ^1 + 1 bj p . For the neutral 

part we get 1 because af = 0. We put b% +1 = b$ ■ (l + j)^ &g, b^ p+1 = b%. 
The stable part is transformed accordingly. 

Part II. The term which is the sum of quadratic functions of all S^'. We expand it using 
the basis of our functions 3>j"\ and the functions from the stable part of the 
spectrum. The result is included in (3^ p \ f3^ p and the stable function $p S ^(Y, s). 

Part III. All remainders which arise because the formulas for finite p are different from the 
limiting formulas. These remainders were estimated in §6. The result is written as 
a series w.r.t. our basis and the corresponding terms are included in (3j U p \ (3^ p and 
the stable part of the spectrum. 



Finally we have 

(«) 

&Sli = &?2 f i + -Y' , &S B 1, = &S B) 



i,p+i j,p y pj ' i'P +1 j'P 

and the formulas for (3^ p \ (3j? p and &p St \Y, s). This works for p < p m +i- If P — Pm+i, then 
we introduce new variables (rescaling!) 

/>) = h (u) I 1 A —I J 4- /? (u) 

i,Pm + l i,Pm + l-l I X T p m+1 J H 3,Pm + V 

b^ = &<?> . + /5 ( r } • 

J ,Pm+l J ,Pm+l-l O ,Pm+l 

It is our other inductive assumption that 

-pT < $L < a>™ < *SL < ^ 

where < pi < 1 but p x is sufficiently close to I. 

Let = [-Pr +1 >Pr +1 ] and A^ = {(&£>^J = < #L +1 < 

It follows easily from the estimates of (3^ 1? /? ( 7l that a£T +1) C A^ } . If Af, m) = 
tffl) : (b^,bp^ m ) G A^j, then A^ is a decreasing sequence of closed sets. The 
intersection f] m gives us the values of parameters for which 5^ — > oo as p — > oo. 

We make also some shortening of the time interval S^" 1 ). In the formulas for 5^ there 
are several remainders which appear when we replace in all expressions s' and s" by s. We 
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estimate these remainders using the fact that our functions satisfy the Lipschitz condition 
and the Lipschitz constants and the maxima of their values decay as some power of p. We 
choose the interval 5'( m + 1 ) c so that when we consider s G S^ m+1 ^ these remainders do 
not violate the basic inclusion A^ +1 - > C • It is easy to see 5 , ( m+1 ) can be chosen so that 
g(m) y g(m+i) COIls i s ts of two intervals whose lengths decay exponentially. Therefore (~) m 
is an interval of positive length. 

The transformation (b^fl , ) — > [b^ 1 ] , b^\ ) is given by smooth functions and is 

sufficiently close to the identity map. The last step in the renormalization is the replacement 
in all 5^ r \ r < p m +i the variables b^ m , b^ pm by their expressions through &^ m+1 , ^> m+ i- 
The form of 5^ in new variables is the same as before. 

The Choice of Constants 

The main constants which are used in the construction are the following: 

1. A;(°) which determines the position of the domain where v(k, 0) is concentrated. 

2. D 1 is the constant which determines the size of the neighborhood where v(k, 0) is 
concentrated. 

3. pi determines the size of the neighborhood where the main parameters b^\ b^ vary. 

4. D2 is the constant which determines the possible size of perturbations in the form 
of v(k,0). 

5. Ai is the power which gives the estimation of the decay of h r in the domain B. 

6. A2 is the parameter which determines the size of the first part of the procedure. 

7. e determines the values of p where we make the renormalization. 

The value of should be sufficiently large. All estimate of the remainders which appear 
during the first half of the procedure are less than const 1 , They should be so small that the 

estimates of all pffi, Pj^ are muc h smaller than p\. On the other hand, p\ should be small 
but not too small. It should be small in order to make the quadratic part of our formulas 
smaller than the linear part. However p cannot be too small in order that we could choose 



52 



the next interval [— p m+1 , p m+l \. This can be achieved by the choice of k^°\ The parameter 
A2 should be small. In this case the estimates of all corrections are easier. However, after A2 
is chosen the value of k^ can be taken sufficiently large depending on A2. The parameter 
Ai can be arbitrarily large in order to make the perturbation arbitrarily small. The value of 

.Di determines the estimates in the domain B which decay as -, — W. We choose D\ so that 

J (fc(°)) 1 

Ai > \. The value of e is chosen small so that we can write with a good precision the action 
of the linearized renormalization group. 



§10. Critical Value of Parameters and Behavior 
of Solutions near the Singularity Point 

We return back to the first formulas: 

t 

v A (k,t) = exp{-t\k\ 2 }A-v(k,0) + / exp {-(* - s)|A;| 2 } • ^ A p h p (k, s) ds 

p>1 

or 



va 



t 

(k,t) = exp{-t\k\ 2 }A- v{k,0) + / exp {- (t - s)\k\ 2 } ■ ^2A p g p {k^,s)ds. (41) 

n P>1 



Our construction gives us the interval S = f] n on the time axis such that for each 
t G S we can find the values of parameters bj = bj(t), 1 < j < 4 and 6^ = bp\t), 
1 < f < 6 such that we have the representation (31) with 8^ — > as r — > 00. It is easy to 
see that A cr (t) = A _1 (t). If so then A p h p (k,t) is concentrated in the domain with the center 
at ^7=^ having the size 0(y/p) and there it takes values 0(p). This immediately implies that 
at t the energy is infinite. 

Consider t! < t. It is important to investigate the behavior of E(t') and the enstrophy 
Q(t') of the same solution with A = A cr (t) when t' is close to t. Denote At = t — t'. It 
follows easily from the proof of the main result that A(t')/ &m = (l — CiAt + 0(At)) for some 
constant C v Since A p cr ■ (A(t')) p = A p r -(A(t)) p - (A(t')/ m ) p = (1 - C At + o(At)) p . It is clear 

"j^ J^Q^ "i~V> o f Drm c in /^^^T\ o vo r»1 acd i~r\ 00 r»Vi r\4- Vi or ~fv w nrt <"" /O I At) j TTVm'* nr\ \ \, ^( At) 



the terms in (|4*T|) are close to each other for p < O ^ &l — )■ For p >> ^J. — the 
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product AP. r (A(t')) p tends exponentially to zero and dominates other terms of the expansion. 
Therefore it is enough to consider \k\ < O — j and in this domain the solution grows as 

3 1 

| k\ 2 . The factor \k\ 2 appears because for any k the values for which the terms in (j41j) give the 
essential contribution to the solution belonging to an interval of the size = 0(*Jp). 



From this argument it follows easily that E(t') ~ a? — j and Q(t') ~ 
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Appendix I. Hermite Polynomials and their basic properties 



Take a > and write 

He^\x) = (-l) n e^-^e-^, n > 0. 

It is clear that He^ (x) = a n x n + • • • , where dots mean terms of smaller degree. We 
shall call He^ the n-th Hermite polynomial. It is clear that He^\x) = 1, Hef\x) = 
ax, He^\x) = a 2 x 2 — a and so on. In general, Hen\x) = a% Hen\\f(Jx). It is easy to 
check that 

axHe^(x) = He { ° ] +1 {x) + anHe^x). (42) 

The Fourier transform of He m ( x ) e ~^~ is (i\) m e~2\ . This implies the formula for 
convolution: 

J H<&l(x - yy-^^-Het\{y)e-^^dy = He% +m2 (x) e -^ ^ (43) 



/ \ ax 

Take positive 71, 72, 71+72 = 1 and consider the convolution of He y m [{-j=)e 271 



j / \ ox 1 A 72 

\ 2^1 an< ^ ^em 2 (^)e • w^j- Their Fourier transforms are (iA^/Ti)™^ -- ^ - and 

A 2 72 Hil m 2 

(iA^/72) m2 e~^~ respectively. The product of these two functions is 7^ 7 2 2 

m l m 2 , s CTX 2 

Therefore the convolution is 7^ • 7 2 2 He m [ +m2 {x)e a - . 



Hil !!i2. \2 
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